Convergence behavior of variational perturbation expansion — 
A method for locating Bender- Wu singularities 
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Variational perturbation expansions have recently been used to calculate directly the strong- 
coupling expansion coefficients of the anharmonic oscillator. The convergence is exponentially fast 
with superimposed oscillations, as recently observed empirically by the authors. In this note, the 
observed behavior is explained and used to determine accurately the magnitude and phase of the 
leading Bender-Wu singularity which is responsible for the finite convergence radius in the complex 
coupling constant plane. 



Variational perturbation theory yields uniformly and exponentially fast converging expansions for many quantum 
mechanical systems The convergence was proved for the anharmonic integral and the quantum mechanical 

anharmonic oscillator in several papers for a finite coupling strength g of the anharmonic term Recently, this 
theory has been adapted to calculating the coefficients of strong-coupling expansions || . For the ground-state energy 
of the anharmonic oscillator, the first 23 coefficients of the series 
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were calculated with 20 digits precision. Thereby an interesting phenomenon was observed: the approach to the 
asymptotic values of a n . is exponentially fast with periodic modulations, as shown in Figs. ^ and ^| calling for a 
revision of the above-quoted proofs which do not apply to the strong-coupling limit g — > oo. The purpose of this 
note is to explain the observed behavior theoretically and to show that it leads to a rather precise estimate of the 
largest Bender-Wu singularities in the complex-g-plane which determine the convergence radius of the strong-coupling 
expansion. 

The potential of the anharmonic oscillator is 



V(x) 



3 4 



( W 2 , 5 >0). 



The Rayleigh-Schrodinger perturbation theory yields a power-series expansion 
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where E^ are rational numbers 1/2, 3/4, -21/8, 333/16, -30885/128, ... , obtained from the recursion relations of 
Bender and Wu 0. 

As is well known, the series (|^) cannot be used for an evaluation of the energy since it has a zero radius of convergence 
due to the factorial growth of the coefficients . The recently developed variational perturbation theory || converts 
the divergent series (|^) into an exponentially fast convergent one. The procedure goes as follows (see Section 5.13 of 
Refs. [||). First, the harmonic term of the potential is split into a new harmonic term with a trial frequency fl, and 
a remainder: 
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After rewriting 



V(x) = ^-x 2 + 9 -{-2ax 2 /VL + x 4 ), 
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we perform a perturbation expansion in powers of g at a fixed a = fl(fl 2 — uj 2 )/g, 
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E N (g,a) = njr £ l 0) (<r)(^)\ (6) 
1=0 ^ ' 

where g = g/Q 3 is a dimensionless reduced coupling constant. 

The calculation of the new series up to a specific order N requires only little work, being easily obtained from the 
ordinary perturbation series (||) by replacing to by yjQ 2 — gcr/Q, and by reexpanding (||) in powers of g up to the Nth. 
order. This yields the reexpansion coefficients 

£ rw=E<( (1 7_ 3 f /2 )(-4^- (7) 

j=0 v / 

The truncated power series Wjy(g,Q) = E^ (<?,cr) is certainly independent of Q in the limit N — > oo. At any 
finite order, however, it does depend on £1 via a, the approximation having its fastest speed of convergence where it 
depends least on Q. If we denote the order-dependent optimal value of f2 by f2j\r, the quantity fijv) is the new 

approximation to E(g) j7|. Introducing the reduced frequency lu — ui/Q, the approximation can be written as 

W N = {g/g) 1/3 w N {g,Lo 2 ). (8) 



From the approximate energies Wn it is easy to derive simple formulas for the coefficients of the strong-coupling 
expansion. We simply expand the function wjv(g,w 2 ) in powers of (D 2 = (g/oj 3 )~ 2 / 3 g 2 / 3 and find 
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with the coefficients 
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Here tojy (a, 0) denotes the n'th derivative of WN(g, & 2 ) with respect to lu 2 at a) 2 = 0, which can easily be calculated 
using 
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The optimal value of fi^r has the TV-dependence (see Ref. [^J and the first of Refs. 
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where the coefficient c is c = 0.186 047 272 .... With this f2jv, we obtain the exponentially fast approach to the exact 
limit as shown in Figs. |l| and g. The exponential falloff is modulated by oscillations. 

To explain this behavior we recall that the ground state energy satisfies the subtracted dispersion relation 
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where disc E^°\g') denotes the discontinuity across the left-hand cut in the complex <?-plane (below minus above). 
An expansion of the integrand in powers of g yields the perturbation series @ . The reexpanded series @ is obtained 
from (|^) by the above replacement of to — ► 51(1 — eg) 1 / 2 and a reexpansion in powers of g. 

There is a simple way of obtaining the same reexpansion from the dispersion relation (I13) . Introducing the dimen- 
sionless coupling constant g = g/to 3 , above replacement amounts to 
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Since Eq. (jjj 
replacement 



represents an energy, it can be written as us times a dimensionless function E^\g). Apart from the 
14]) in the argument, it receives an overall factor fl/w = (1 — erg) 1 / 2 . If we introduce the reduced energy 
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E(g) = E(g)/Q, which depends only on the reduced coupling constant g, the dispersion relation ( fl3| ) for E^(g) turns 
into 



E^\g) = {l^agf" 
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The resummed perturbation series is obtained from this by an expansion in powers of g/A up to order N . 

Note that only the truncation of the expansion causes a difference between the two expressions (13) and (1151), since 
g and g are the same numbers, as can be verified by inserting g = g/fl 3 and a into the right-hand side of (^4[). 

To find the reexpansion coefficients we observe that the expression ( |l5| ) satisfies a dispersion relation in the complex 
g-plane. If C denotes the cuts in this plane and disccE(g) is the discontinuity across these cuts, the dispersion relation 
reads 
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We have changed the argument of the energy from g to g since this will be the relevant variable in the sequel. 

When expanding the denominator in the integrand in powers of g/A, the expansion coefficients ej ^ are found to be 
moment integrals with respect to the inverse coupling constant 1/g: 
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In the complex g-plane, the integral ( |l5| ) has in principle cuts along the contours Ci,Cj,C2,C2, and C3, as shown 
in Fig. [| The first four cuts are the images of the left-hand cut in the complex g-plane; the curve C3 is due to the 
square root of 1 — og in the mapping (|lj) and the prefactor of (fl5|). 

Let D(g) abbreviate the reduced discontinuity in the original dispersion relation (filf): 



D{g) = discE^ia) = 2ilmE {0) (3 - ir)), d < 0- 
Then the discontinuities across the various cuts are 

disc E^{g) = (1 - ag) l / 2 D(g(l - a.9)- 3 / 2 ), 



disc E^\g) = -2i{ug- 1) 1/2 
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For small negative g, the discontinuity is given by the semiclassical limit (see 



Chapter 17): 



D{g) 



-2i 



=4/3g 



7T V -3g 



(18) 
(19) 

(20) 
(21) 



We denote by ejj.° (Cj) the contributions of the different cuts to the integral ( |l7|) for the coefficients. After inserting 
( pl| ) into Eq. (|l^), we obtain from the cut along C\ the semiclassical approximation 



ef{C l )^~2A k 
For the fcth term Sk of the series this yields an estimate 
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where fk{l) is the function of 7 = og 

fkil) 
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For large k, the integral may be evaluated via a saddle point approximation. At the extremum, 7 — > 7^ = — 4er/3fc, 

k — >oo 

fk{l) has the value 

f k — » fclog(3fc/4e ( 7)-2a. (25) 

— >oo 

The constant — 2a in this limiting expression arises when expanding the second term of Eq. (|24|) into a Taylor series, 
(47/37) (1 — 7) 3 ^ 2 — 4cr/37fe — 2cr + . . . . Only the first two terms survive the large-fc limit. 
Thus, to leading order in k, the fcth term of the reexpanded series becomes 

S k oce-^(^) k (^) k . (26) 



The corresponding reexpansion coefficients are 



4 0) °ce- 2 ^f. (27) 



They have the remarkable property of growing in precisely the same manner with k as the initial expansion coefficients 
except for an overall suppression factor e _2<T . This property was discussed in Ref. |2). 
In order to estimate the convergence of the variational perturbation expansion, we note that ag = 1 — 1/fl 2 . For 
large O, this expression is smaller than unity. Hence the powers (crg) k alone yield a convergent series. An optimal 
reexpansion of the energy can be achieved by choosing, for a given large maximal order N of the expansion, a 
parameter a proportional to N: 

a w <r N = cN. (28) 
Inserting this into (|24|), we obtain for large k = N 

Ac 
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In{i) « N 
The extremum of this function lies at 

1+^(1_7)V2 (1+ 1 7)=0 . (30) 

The constant c is now chosen in such a way that the large exponent proportional to N in the exponential function 
e /iv(7) cl U e to the first term in (^9|) is canceled by an equally large contribution from the second term, i.e., we require 
at the extremum 

M7) = 0. (31) 

The two equations ([50|) and ( |3l| ) are solved by 

7 = -0.242 964 029 973 520 .. ., c = 0.186 047 272 987 975 ... . (32) 

In contrast to the extremal 7 of Eq. ( p4| ) which dominates the large-fc limit, the extremal 7 of the present limit, in 
which k is also large but of the order of N, remains finite (the previous estimate holds for k 3> N). Accordingly, the 
second term (43/37) (1 — 7) 3 ^ 2 m /jv(7) contributes in full, not merely via the first two Taylor expansion terms of 
(1 -7) 3/2 , as it did in (123). 

Since /a? (7) vanishes at the extremum, the A^th term in the reexpansion has the order of magnitude 

S N (x (cr N g N ) N = ^1 - . (33) 

According to (|2|), the frequency flpf grows for large N like 

tlN~atf l g l '*~{eNg)W. (34) 
As a consequence, the last term of the series decreases for large N like 
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This estimate does not yet explain the convergence of the variational perturbation expansion in the strong-coupling 
limit observed in Figs. ^ and ||. For the contribution of the cut C\ to Sn, the derivation of such a behavior requires 
including a little more information into the estimate. This information is supplied by the empirically observed property, 
that the best f^-values lie for finite N on a curve (see Chapter 5 in 
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Thus the asymptotic behavior (|28| ) receives, at a finite TV, a rather large correction. By inserting this <jn into /jv(t) 
of (129), we find an extra exponential factor 
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This reduces the size of the last term due to the cut SW(Ci) in ( p5[ ) to 

5 JV (C 1 )oce-[ 9 - 7 +^)" /3 ] Jvl/3 ) 



(37) 
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which agrees with the convergence seen in Figs. and |. 

Note that there is no need to evaluate the effect of the shift in the extremal value of 7 caused by the correction 
term in (|36|), since this would be of second order in I/N 2 / 3 . 

How about the contributions of the other cuts? For Cj, the integrals in ( JT7| ) runs from g = — 2/ a to —00 and 
decrease like {—2/a)~ k . The associated last term SV(Ci) is of the negligible order e _A " ogW . For the cuts C 2 ,2,3, the 
integrals ( |l7j ) start at g = 1 ja and have therefore the leading behavior 
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This implies a contribution to the iVth term in the reexpansion of the order of 



(39) 



(40) 



which decreases merely like ( pq ) and does not explain the empirically observed convergence in the strong-coupling 
limit. As before, an additional information produces a better estimate. The cuts in Fig. [| do not really reach the 
point ag = 1. There exists a small circle of radius Ag > in which E^°\g) has no singularities at all. This is a 
consequence of the fact unused up to this point that the strong-coupling expansion (0) converges for g > g s . For the 
reduced energy, this expansion reads: 
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The convergence of (|lj) for g > g s implies that ( |4l|) converges for all ag in a neighborhood of the point ag = 1 with a 
radius 



A(ag) ~ ( -4- 
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where g s = g s /uj 3 . For large A^,A(crg) goes to zero like l/(N\gs\c) 2 / 3 . Thus the integration contours of the moment 

integrals ( |l7j ) for the contributions ef\Ci) of the other cuts do not begin at the point ag = 1, but a little distance 
A(ag) away from it. This generates an additional suppression factor 
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Let us set — g s = \g s \ exp(iip s ) and x s = (—g/g s ) 2 ^ — — \x s \ exp(i#), and introduce the parameter a = 1 / [ | <7 S | c] 2 / 3 . 
Since there are two complex conjugate contributions we obtain, for large N a last term of the order of 

Sn(C 2X3 ) » e~ Nl/3acose cosiN^asmO). (44) 

By choosing 

|&|~ 0.160, 6>- -0.467, (45) 

we obtain the curves shown in Figs. || and [|, which agree very well with the observed Figs, [j] and [|. Their envelope 
has the asymptotic falloff e - 9 - 23Arl/3 . 

Let us see how the positions of the largest Bender- Wu singularities compare with what we can extract directly 
from the strong-coupling series (Q) up to order 22 Q. For a pair of square root singularities at x s = — \x s \ exp(±i0), 
the coefficients of a power series ^a n x n have the asymptotic ratios R n = a n+ i/a n — = — cos[(n + 1)9 + 
5}/\x s \ cos(n9 + S). Plotting the ratios R n for the coefficients a n in Fig. || we see that for large n, they are well 
reproduced by R™ if we choose 

|z 8 | = 1/0.117, 9 = -0.467, (46) 

with an irrelevant phase angle S = —0.15. The angle 9 is in excellent agreement with the value found in (ff5|). From 
\x s \ we find \g s \ = 4|l/x s | 3 / 2 = 0.160, again in excellent agreement with (|45|). 

Note that this convergence radius is compatible with the heuristic convergence of the strong-coupling series up to 
order 22, as can be seen in Fig. ^ by comparing the series with the exact curve. 



While thi s work was in progress, we received a Genova preprint by R. Guida, K. Konishi, and H. Suzuki ( hep- 1 
th/! j505084 ) in which the exponentially fast convergence found in Ref. [|| is proved rigorously by reducing it to the 



known convergence of the strong-coupling expansion. 
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FIG. 1. Exponentially fast convergence of the Nth approximants for a to the exact values. The plotted quantity is 
A N = |(ao)jv) - a \. 




FIG. 2. Oscillatory behavior around the asymptotic approach of Qo to its exact value as a function of the order N of the 
approximant (open circles are for odd N, filled circles for even N). 



FIG. 3. Cuts in the complex g-plane whose moments with respect to the inverse coupling constant determine the reexpansion 
coefficients. The cuts inside the shaded circle happen to be absent due to the convergence of the strong-coupling expansion for 
g > g s (from Ref . §) . 
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FIG. 4. Theoretically obtained convergence behavior of the Nth approximants Sn for ao, to be compared with the empirically 
found behavior in Fig. 1. 
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FIG. 5. Theoretically obtained oscillatory behavior around the asymptotic approach of a to its exact value as a function 
of the order TV of the approximant, to be compared with the empirically found behavior in Fig. 2, averaged between even and 
odd orders. 
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FIG. 6. Comparison of ratios R n between successive expansion coefficients of the strong-coupling expansion (dots) with the 
ratios of the pynansinn nf a «nr>erposition of two singularities at g = 0.156 x exp(±0.69) (crosses). 
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FIG. 7. Strong- Coupling Expansion for the ground-state energy in comparison with the exact values and the perturbative 
results of 2nd and 3rd order. The convergence radius in 1/g is larger than 1/0.2. 
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